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We consider a particle in the harmonic approximation coupled linearly to an environment, mod- 
eled by an infinite set of harmonic oscillators. The system (particle-environment) is considered in 
a cavity at thermal equilibrium. We employ the recently introduced notion of renormalized coordi- 
nates to investigate the time evolution of the particle occupation number. For comparison we first 
present this study in bare coordinates. For a long ellapsed time, in both approaches, the occupation 
' number of the particle becomes independent of its initial value. The value of ocupation number of 

C^' the particle is the physically expected one at the given temperature. So we have a Markovian pro- 

cess, describing the particle thermahzation with the environment. With renormalized coordinates 
' no renormalization procedure is required, leading directly to a finite result. 

^. . . 

A thermahzation process occurs in some cases for a system of material particles coupled to an environment, in the 
sense that after an infinitely long time, the matter particles loose the memory of their initial states. This study is, in 
general, not easy from a theoretical point of view, due to the complex non-linear character of the interactions between 
■ the matter particles and the environment. To get over these difficulties, linearized models have been adopted. An 
' account on the subject of the evolution of quantum systems on general grounds can be found in [l|, 0, H, @]- 
. Besides, the main analytical method used to treat these systems at zero or finite temperature is, except for a few special 
' cases, perturbation theory. In this framework, the perturbative approach is carried out by means of the introduction 
, of bare, non-interacting objects (fields, to which are associated bare quanta), the interaction being introduced order 
by order in powers of the coupling constant. 

In spite of the remarkable achievements of the perturbative methods, however, there are situations where they 
cannot be employed, or are of little use. These cases have led to attempts to improve non-perturbative analytical 
methods, in particular, where strong effective couplings are involved. Among these trials there are methods that 
^ , perform resummations of perturbative series, even if they are divergent, which amounts in some cases to extending 
' the weak-couplin g reg ime to a strong-coupling domain. One of these methods is the Borel resummation of perturbative 

rs ; series 0, i i [isnnriii • 

• In this paper we follow a different non-perturbative approach: we investigate a simplified linear version of a particle- 
field or particle-environment system, where the particle, taken in the harmonic approximation, is coupled to the 
reservoir, modeled by independent harmonic oscillators @. We will eniploy, in particular, dressed states and 

renormalized coordinates, introduced in [T3] and already employed in [Tsl . [igI. IITI. [isj . Using this method non- 
perturbative treatments can be considered for both weak and strong couplings. A linear model permits a better 
understanding of the need for non-perturbative analytical treatments of coupled systems, which is the basic problem 
underlying the idea of a dressed quantum mechanical system. Of course, the use of such an approach to a realistic 
non-linear system is an extremely hard task, while the linear model provides a good compromise between physical 
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reality and mathematical reliability. The whole system is supposed to reside inside a spherical cavity of radius R in 
thermal equilibrium at temperature T = I3~^ . In other words, we consider the spatially regularized theory (finite 
K) at finite temperature. The free space case is obtained by suppressing the regulator, (i? oo). For a detailed 
comparison between this procedure and the one considering an a priori unbounded space, see [l4i] . 



I. THE MODEL 



Let us start by considering a particle approximated by a harmonic oscillator, having bare frequency wq, linearly 
coupled to a set of N other harmonic oscillators, with frequencies ojk, k = 1,2, . . . , N . The Hamiltonian for such a 
system is written in the form, 



"4 



N 

fe=i 



{pI 



2 2\ 



N 

fc=i 



(1) 



leading to the following equations of motion. 



N 



qi + qi 



= ^c^q^{t) 

i=l 

= c,qo{t). 



(2) 
(3) 



In the limit ^ oo, we recover our case of the particle coupled to the environment, after redefining divergent 
quantities, in a manner analogous to mass renormalization in field theories. A Hamiltonian of the type ([T]) has been 
largely used in the literature, in particular to study the quantum Brownian motion with the path-integral formalisni 
^ . It has also been employed to investigate the linear coupling of a particle to the scalar potential [Tj, [3, [Hj [13) fig - 
The Hamiltonian ([T]) is transformed to principal axis by means of a point transformation. 



N 



N 



= ^^l^Qr, Pl^=^tl,Pr; 

r=0 r=0 

fi = (0,{fc}), A; = 1,2,. ..,7V; r = 0,...N, 



(4) 



performed by an orthonormal matrix T — {f^). The subscripts /i = and fi = k refer respectively to the particle and 
the harmonic modes of the reservoir and r refers to the normal modes. In terms of normal momenta and coordinates, 
the transformed Hamiltonian reads 



1 ^ 



(5) 



where the fi^'s are the normal frequencies corresponding to the collective stable oscillation modes of the coupled system. 
Using the coordinate transformation Q in the equations of motion and explicitly making use of the normalization of 



the matrix (r), Eu=o(^m)^ = 1' get 



Cfe 



r 



with the condition 



N 

E 

k=l 



N 



k=l 



(6) 



(7) 



We take = 77(0;^)", where ry is a constant independent of k. In this case the environment is classified according to 
M>l,u=l, oru<l, respectively as supraohmic, ohmic or subohmic. This terminology has been used in studies of 
the quantum Brownian motion and of dissipative systems For a subohmic environment the sum in Eq. ^ 

is convergent in the limit N ^ 00 and the frequency luq is well defined. For ohmic and supraohmic environments, this 
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sum diverges for N oo. This makes the equation meaningless, unless a renormalization procedure is implemented. 
From now on we restrict ourselves to an ohmic system. In this case, Eq. ([7]) is written in the form 



u^l-s.'~nl^^-nlY^-^^, (8) 

k=l '^k ^^r 

where we have defined the counterterm 

Slu^ = Nri^. (9) 

There are + 1 solutions of fi^, corresponding to the + 1 normal collective modes. Let us for a moment suppress 
the index r of fi^. If Wq > (5w^, all possible solutions for are positive, physically meaning that the system oscillates 
harmonically in all its modes. If ujq < SuP', then a single negative solution exists. In order to prove this let us define 
the function 

^ 1 

= ul -V?- n^V? Y: -2^^, (10) 

fe=i '^k " 

so that Eq. dH]) becomes / (fi^) = 0. We find that 

/(ri'^)^oo as — oo and /(O) = cjq — Suj'^ < 0, 

in the interval (— oo, 0]. As / {^^) is a monotonically decreasing function in this interval, we conclude that / (f^^) = 
has a single negative solution in this case. This means that there is a mode whose amplitude grows or decays 
exponentially, so that no stationary configuration is allowed. Nevertheless, it should be remarked that in a different 
context, it is precisely this runaway solution that is related to the existence of a bound state in the Lee-Friedrichs 
model. This solution is considered in the framework of a model to describe qualitatively the existence of bound states 
in particle physics (19j . 

Considering the situation where all normal modes are harmonic, which corresponds to the first case above, Wq > Suj'^, 
we define the renormalized frequency 

uj'^ =ujI- Suj'^ = lim {uol - Nrj^), (11) 

N^OO 

in terms of which Eq. ^ in the hmit N ^ oo becomes, 

k=l « 

In this limit, the above procedure is exactly the analog of the mass renormalization in quantum field theory: the 
addition of a counterterm —Suj'^qQ allows one to compensate the infinity of luq in such a way as to leave a finite, 
physically meaninful renormalized frequency uj. This simple renormalization scheme has been introduced earlier [20| . 
Unless explicitly stated, the limit — > oo is understood in the following. 
Let us define a constant g, with dimension of frequency, by 



where Acu ~ nc/R. The environment frequencies ujk ^"^^ given by, 



7TC 

'^k = k—, fc = l,2,..., (14) 



where R is the radius of the cavity that contains the whole system. Then, using the identity 

(15) 



1 1 

k=l 



1 vr 

cot (ttu) 



Eq. p2|) can be written in a closed form: 



The solutions of the above equation with respect to f2 give the spectrum of eigenfrequencies fir corresponding to the 
collective normal modes. 

In terms of the physically meaningful quantities fi^ and Cu, the transformation matrix elements turning the particle- 
field system to the principal axis are obtained. They are 



f"^ — 
— 



(17) 



These matrix elements play a central role in the quantities describing the system. 



II. THE THERMALIZATION PROCESS IN BARE COORDINATES 

We now consider the thermalization problem using bare coordinates. For the model described by Eq. ([1]) this 
problem was addressed in an alternative way in [21j with the canonical Liouville-von Neumann formalism. We 
consider the initial state described by the density operator, 



p{t = 0) = po ® P/3 



(18) 



where po is the density operator of the particle, that in principle can be in a pure or in a mixed state and p/j is the 
density operator of the thermal bath, at a temperature that is, 



P/9 = ^ exp 



(19) 



with Zp = rifcLi being the partition function of the reservoir, and 



Trfc 



-l3LOk{alak+l/2) 



The creation and annihilation operators given by 



1 



2 sinh 



(^) 



(20) 



(21) 
(22) 



where cZj^ — {uj,ujk). The thermalization problem is addressed by investigating the time evolution of the state p(t). 

The thermalization problem concerns the time evolution of the initial state to thermal equilibrium. The subsystem 
corresponding to the particle oscillator is described by an arbitrary density operator po- As we will show, the 
expectation value of the number operator corresponding to particles will evolve in time to a value that is independent 
of the initial density operator po, the dependence will be exclusively on the mixed density operator corresponding to 
the thermal bath. 

Our aim is to obtain expressions for the time evolution of the expectation values for the occupation number and 
in particular for the one corresponding to particles. We will solve the problem in the framework of the Heisenberg 
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picture. It is to be understood that when a quantity appears without the time argument it means that such quantity 
is evaluated aA, t = 0. The Heisenberg equation of motion for the annihilation operator a^(t) is given by 



Due to the linear character of our problem, this equation is solved by writing a^(t) as 

oo 



(23) 



(24) 



where all the time dependence is in the c -rawm&er functions B^^{t). Then, Eq. (|23p reduces to the following coupled 
equations for B^„{t): 



fc=i 



Bf,k{t) + ojlB^kit) - B^n{t) V Tjujk = 



fe=i 



(25) 



(26) 



These equations are formally identical to the classical equations of motion, Eqs. ([2]) and ([3|), for the bare coordinates 
Qfj,. Then we decouple Eqs. and (pS)) with the same matrix {i^} that diagonalizes the Hamiltonian Eq. ([1]). In 
an analogous manner, we write B^^{t) as 



r=0 



such that from Eqs. (|25p and (pS)) . we obtain the following equations for the normal-axis functions CJ^{t), 
which gives the solution 

Then substituting this expression into Eq. (P?)) we find 



(27) 



(28) 



(29) 



r=0 



The time independent coefficients a]^, 6^ are determined by the initial conditions at t = for Bfj_^{t) and B^,y{t). 
From Eqs. (PTjl and (IM)) we find that these initial conditions are given by 



B., 



(30) 



Using these equations, we obtain for and 



(31) 
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We write a^(i) and aj^(i) in terms of and using Eqs. (I2ip . ((22)) and p4l) . 

oo 

= ^ (aAii/(t)ai/ + /3pi.(t)a|;) , (33) 

1^=0 

oo 

= E(/3^-W«- + "M^W^t) , (34) 
where afi^{t) and P^^it) are the Bogohubov coefficients given by, 

a^,{t) = -^B^„(t) - i, r±B^,{t) (35) 



and 



/3^,(i) = -^i3^^[t) + t.^B^^it) . (36) 

' ZLsJ,j V 



Using the definition of B^iy{t) we get 

+ [{ftr ~ Lu,,)e'''^* + {ilr + w^)e-*"'-*] } , (37) 



and 

r—O V V 

- [{Qr - u;f,)e'''^' + {a,. + u;^)e-''''-'] } . (38) 
Now we study the time evolution of n^(t), the expectation value of the number operator -/V^(t) = (t)a^(<), that is, 

nf,{t)^Ti[al{t)a^{t)po(E)pp] . (39) 
Using the basis |no, ni, 712, ...un) we obtain. 



= E + \P,Am + E l/^M. Wl' ' (40) 

where 



u=0 iy=0 



00 



no = 2_^n{n\po\n) (41) 

n=0 

is the expectation value of the number operator corresponding to the particle and the set {uk} stands for the thermal 
expectation values corresponding to the thermal bath oscillators, given by the Bose-Einstein distribution, 

nk = -TT-^ • (42) 

In Eq. ([40l) there appears a term that does not depend on the temperature of the thermal bath. This term has its 
origin in the instability of the initial bare vacuum state, |0, 0, 0). To see this, we compute the expectation value of 
the time dependent number operator iV^(t) = aj((t)a^(t) in this vacuum state. Thus all the terms containing operators 
different from the identity give a zero contribution. The only term, that gives a non-zero contribution comes from the 
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normal ordering and is just the last one in Eq. (|40)) . This term leads to the creation of excited states (particles, in a 
field theoretical language) from the initial unstable bare vacuum state. 

We are interested in evaluating the expectation value of the number operator corresponding to particles. Thus 
taking /i = in Eq. and using Eq. ([32]), we obtain 



no{t) = [Mt)\^ + \l3oo{t)\^]no + J2[\aokit)f + \Mt)\' 



k=l 



k=l 

where the coefficients of this expression are [2l| , 



and 



such that 



"00 (0 
/3oo(i) 
aokit) 



IQuJK 

■Kgt/2 



(2w + 2k — ii:g) e 



2 iKt 



nge 



8UJK 



- {2lu -2k- iirg) e 
{irg - 2iK) e'"*] 



2a; (w^ — uj^ + ingujk) 
{2n + 2u}- ing) 
(2k — 2ujk — iirg) 



+ 



LUk \/2gAu 



u) 4k 
{2lIi -2k- ing) 
(2k + 2ujk + ing) 



-7rgt/2 



Mt) 



ujk (wfc - uj) y/gAuj e"^'' * juj^ ^/2gKuj 



2uj (uj'j, — up- — inguJk) 
{2il! + 2k- ing) .^^ 



u) Ak 



(2k + 2ujk — ing) 



{2uj — 2k — ing) 

-f 

(2k — 2ujk + iT^g) 



(43) 

(44) 
(45) 

(46) 



(47) 



K = v/'^2_^2^2/4_ (48) 

The parameter k measures the intensity of the interaction: if k^ >> 0, i.e g << 2uj/n, we are in the weak coupling 
regime. On the contrary if k^ << 0, i.e g >> 2uj/n, the system is in the strong coupling regime. Here we will restrict 
ourselves to the weak coupling regime. This case includes the important class of electromagnetic interactions, g = acu, 
with a being the fine structure constant a = 1/137 p^ . 

In the continuum limit Auj —^ 0, sums over k become integrations over a continuous variable uj and we obtain for 
no(t). 



where 



no{t) 



-Trgt 



2 

+ — ^ {2u}^ - n'^g^) cos(2«-i) 



n^g^K 

— - — sin(2Ki) 



no 



n^g2e—9t 



_„ „ \2lu^ + (2uj^ - n^g^) cos(2 Kt) -2ngK sm( 2 Kt)] 



duj 



F{uj,uj,g,t) 
- 1) 



+ Gi^,^,9,t) 



(49) 



F{u),uj,g,t) 



w(tj2 + w^) 



uj''^ — uj'^)'^ + n'^g'^oj 



2„2,,,21 



-Trgt 



4^2 



-[Au)^ - n^g^ cos(2«-i) 



'2ngK 



(w2 + (1,2) 

4070;^ 



sm{2Kt)\ 



,-7rgt/2 



- [2k cos{uot) cos(Kt) 



„ „, sm(ujt) sin(Kt) —nq—^ ^ cosfwt) sm(Kt) 

2 4. -,2^ ^ ^ ' (w2 + (2)2) ^ ' ^ 



(cj2 4- (1,2) 



(50) 
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and 



--^ 9 tt/ cos(2Ki) - 27r.gK-^ rr sin(2Kt) 

(w — w)"^ (cj — Uj) 

p-7rgt/2 

- [2k cos(a;t) cos(Kt) — 2Li; sin([x)f;) sin(Kt) 



K 

— TTff 1^ j cos(a;t) sin(Kt)] . (51) 

It is to be noticed that the second and the third hnes in Eq. ((49|l are independent of the initial distribution. Also the 
integral over G{u!,lu, g,t) is logarithmically divergent. We can understand the origin of these terms in the following 
way: suppose that initially, in the absence of the linear interaction, we prepare the system in its ground state, that is, 
at i = we have |0, 0, 0). Then, we can compute, in the Heisenberg picture, the time evolution for the expectation 
value of the number operator corresponding to the particle, that is (0,0, 0|aQ(i)ao(t)|0, 0, ...,0). We obtain, 

oo 

(0,0,...,0|aS(i)ao(t)|0,0,...,0) = |/3oo(OI' + El^ofcWI' ' (^2) 

fe=i 

which in the continuum limit gives the second and third lines of Eq. (|49[) . Then, these terms appearing in Eq. (|49|) 
are interpreted as the excitations produced from the unstable bare (vacuum) ground state, as a response to the onset 
of the linear interaction. 

The result above is compatible with some results in Q in the context of quantum dissipative phenomena. In 
this quoted paper, in the zero-temperature situation, the system is represented by a set of harmonic oscillators. A 
detailed justification for representing the environment by a set of harmonic oscillators is given in the appendix C of 
this reference. 

The divergent integral in G{ll!^lo, g,t) can be dealt with by a renormalization procedure. The suppression of this 
term is analogous to the standard Wick-ordering in field theory. Thus we write the following renormalized expectation 
value for the particle number operator. 



- m T^/- ^\ , 9 , F{w,w,g,t) 

no{t) ^ K{u;,g,t) + - duj ^^^^ _ (53) 



where 



Kiio,g,t) 



+ ^^-^ (2w2 - TT^g^) cos(2Kt) - ^ ^ ^ sin(2Kt) 



no 



16a)2«;2 



[2lo^ + (2w^ - TT^g^) cos(2Kt) - 27r5Ksin(2Kt)] 



(54) 



In the limit t — > oo, no{t) has a well defined value, that is, the system reaches a final equilibrium state. Also, 
since K{u), g,t ^ oo) — > 0, this final equilibrium state is independent of tiq. The equilibrium expectation value of the 
number operator corresponding to the particle is independent of its initial value, and the only dependence is on the 
initial distribution of the thermal bath, that is, the particle thermalizes with the environment. Before the interaction 
enters into play for i < 0, n{t < 0) = rto, then we have that K{uj,il},g,t < 0) = 1. Taking i = in Eq. (|5T 



we obtain K{uj,uj,g,t — 0) 



T a 



no 



TV g 



[2u}'^ + (2cj2 - TT^g^)] . Thus K{u), g, t) is a 



discontinuous function of t] the discontinuity appearing just at i = 0. From the physical standpoint this discontinuity 
can be viewed as a response to the sudden onset of the interaction between particles and the environment. 

It should be mentioned that a very similar problem from the mathematical point of view, has been studied in [l^. 
In this work the authors study the damped harmonic oscillator under the optics of a a dissipation problem. They 
apply a method that diagonalizes the Hamiltonian of the system and derive the conditions of validity of the rotating 
wave approximation. 
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FIG. 1: Time behavior for no{t) given by Eq. ((53l) for {t > 1), no ^ 1, lj ^ 1, P = 2 and g = 0.1 



Although the integral in Eq. (1531) can not be computed analytically, we can perform numerical calculations, for 
example in Fig. [l]we display the time behavior for no = 1, w = 1, /3 = 2 and g = 0.1; (t > 1). In the next section 
we develop an alternative approach based on the notion of dressed particles. We will find that, in this new realm, no 
renormalization is needed. 



III. DRESSED COORDINATES AND DRESSED STATES 



Let us start with the eigenstates of our system, |no, ni, n2...), represented by the normalized eigenfunctions in terms 
of the normal coordinates {Qr}i 



onin2... 



Tne ^= = = 



where Hn^ stands for the Tis-th Hermite polynomial and Tq is the normalized vacuum eigenfunction. 



(55) 



(56) 



We introduce dressed or renormalized coordinates gg and {q'j} for, respectively, the dressed particle and the dressed 
field, defined by, 

= X! tliV^Qr, (57) 



valid for arbitrary R and where — {uj, uji}. In terms of dressed coordinates, we define for a fixed instant, t = 0, 
dressed states, \ko, ki, K2...) by means of the complete orthonormal set of functions 




(58) 



where q'^ — {qq, g,-}, = {iD, uJi}. Notice that the ground state Fq in the above equation is the same as in Ea. (|55p . 
The invariance of the ground state is due to our definition of dressed coordinates given by Eq. ([57]) . Each function 
V'koki...('Z') describes a state in which the dressed oscillator q'^ is in its K^-th excited state. 
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It is worthwhile to note that our renormalized coordinates are new objects, different from both the bare coordinates, 
q, and the normal coordinates Q. In particular, the renormalized coordinates and dressed states, although both are 
collective objects, should not be confused with the normal coordinates Q, and the eigenstates Eq. ([55|l . While the 
eigenstates (p are stable, the dressed states ip are all unstable, except for the ground state obtained by setting {k^ — 0} 
in Eq. ((58)) . The idea is that the dressed states are physically meaningful states. This can be seen as an analog of 
the wave-function renormalization in quantum field theory, which justifies the denomination of renormalized to the 
new coordinates q'. Thus, the dressed state given by Eq. ([55)) describes the particle in its Ko-th excited level and 
each mode k of the cavity in the Kk — th excited level. It should be noticed that the introduction of the renormalized 
coordinates guarantees the stability of the dressed vacuum state, since by definition it is identical to the ground state 
of the system. The fact that the definition given by Eq. ()57p assures this requirement can be easily seen by replacing 
Eq. (|57)) in Eq. (|55)) . We obtain TqW) Tq{Q), which shows that the dressed vacuum state given by Eq. (|55)) is the 
same ground state of the interacting Hamiltonian given by Eq. 

The necessity of introducing renormalized coordinates can be understood by considering what would happen if we 
write Eq. (l58)) in terms of the bare coordinates 9^. In the absence of interaction, the bare states are stable since 
they are eigenfuntions of the free Hamiltonian. But when we consider the interaction they all become unstable. The 
excited states are unstable, since we know this from experiment. On the other hand, we also know from experiment 
that the particle in its ground state is stable, in contradiction with what our simplified model for the system describes 
in terms of the bare coordinates. So, if we wish to have a nonperturbative approach in terms of our simplified 
model something should be modified in order to remedy this problem. The solution is just the introduction of the 
renormalized coordinates q'^ as the physically meaningful ones. 

In terms of bare coordinates, the dressed coordinates are expressed as 

9m = X!"^''^'^' 



where 



^j2e^t:^,. (60) 



If we consider an arbitrarily large cavity (i? —>■ 00), the dressed coordinates reduce to 



q'o = AQo{uj,g)qo, (61) 
ql = q^, (62) 

with Aoo(w,5) given by. 

In other words, in the limit R 00, the particle is still dressed by the field, while for the field there remain bare 
modes. 

Let us consider a particular dressed state |r5'(0)), represented by the wavefunction ■(/'oo - i(ai)o - (9')- I* describes the 
configuration in which only the dressed oscillator q'^ is in the first excited level. Then the following expression for its 
time evolution is valid fl3 |: 

Kit)) = 5]r''(<)|rr(o)) 

nt) = (64) 

S 

Moreover we find that 

Y^\r^t)f = l. (65) 

V 

Then the coefficients J^^it) are simply interpreted as probability amplitudes. 

In approaching the thermalization process in this framework, we have to write the initial physical state in terms of 
dressed coordinates, or equivalently in terms of dressed annihilation and creation operators and aj| instead of 
and ajj. This means that the initial dressed density operator corresponding to the thermal bath is given by 
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Pl3 = exp 



where we define 



fc=i 



2 



1^, 
2 



(66) 



(67) 
(68) 



Now we analyze the time evolution of dressed coordinates. 



IV. THERMAL BEHAVIOR FOR A CAVITY OF ARBITRARY SIZE WITH DRESSED COORDINATES 

The solution for the time-dependent annihilation and creation dressed operators follows similar steps as for the bare 
operators. The time evolution of the annihilation operator is given by, 



-a;(i) = i [i7,a;(i) 

and a similar equation for a'^it). We solve this equation with the initial condition at i = 0, 



(69) 



which, in terms of bare coordinates, becomes 



N 



r r,r - , ^ 



We assume a solution for a' {t) of the type 



CX3 



j/=0 



Using Eq.((T]) we find. 



r=0 



(70) 



(71) 



(72) 



(73) 



In the present case the time independent coefficients are different from those in the bare coordinate approach, Eq. (|29p . 
The initial conditions for B'^^{t) and iS'^^it) are obtained by setting t = in Eq. ([7^ and comparing with Eq. (|7T|) : 
Then 



OO 

Km - E 

r=0 



tit 



fj. I' 



-t:.t:. 



(74) 
(75) 



Using these initial conditions and the orthonormality of the matrix {f^} we obtain a'^ = 0, 5J[ = if^^J \/2Vlr. Replacing 
these values for a'^ and h'^ in Eq. ([73]) we get 



s;.w = *E 



OO ,r 4-r 



r=0 



(76) 
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We have 

N 



r,i/— 
N 




= E ( + TTTTP'^ ) = E U^i^)^''^ ' (77) 

r,i/— 

where 

oo 

f,At)^Y.^ltle-^''^' . (78) 

For the occupation number n'^{t) = {a'^{t)a'^{t)) we get 

n'^{t) = Tr{a'l{t)a'^{t)p'^®p'p) . (79) 

where p'^ is the density operator for the dressed particle and p'^ is the density operator for the thermal bath, which 
coincides with the corresponding operator for the bare thermal bath if the system is in free space (in the sense of an 
arbitrarily large cavity) [H, . 

To evaluate n'^(t) we choose the basis |no, ni, un) = Il^o \^^^)^ where are the eigenvectors of the number 
operators a^Ja^. From Eq. ([77)1 we get 

oo 
oo 

= Ei^''Wi'«^^«^ + E/;pW^''W«^^«p- (80) 

I/ — l^T^P 

In the basis |n,o, ?^27 ■ ' ') obtain, 

oo 

<(i) = l/Mo(OlX + El/pfe(^)l'<' (81) 

k=l 

where ti'q and nj. are the expectation values of the initial number operators, respectively, for the dressed particle and 
dressed bath modes. We assume that, dressed field modes obey a Bose-Einstein distribution. This can be justified 
by remembering that in the free space limit, R — s- oo, dressed field modes are identical to the bare ones, according 
to Eqs. (j6ip and (j62[) . Now, no term independent of the temperature appears in the thermal bath. This should be 
expected since the dressed vacuum is stable, particle production from the vacuum is not possible. Setting p, — Q va 
Eq. (j8ip we obtain the time evolution for the ocupation number of the particle, 

oo 

"o(i) = l/oo(i)lX+El/o'=(^)lX- (82) 

V. THE LIMIT OF ARBITRARILY LARGE CAVITY: UNBOUNDED SPACE 

In a large cavity (free space) we must compute the quantities /oo(0 fok{t) in the continuum limit to study the 
time evolution of the ocupation number for the particle. Remember that in Eqs. (jl7p . ujk — knc/R, k = 1,2, ... and 
1] = y/2gAuj, with Aw = {uji+i — iLii) — nc/R. When R — > oo, we have Auj — > and AJ7 and then, the sum in 
Eq. ((78|) becomes an integral. To calculate the quantities f^u{t) we first note that, in the continuum limit, Eq. (fT7|) 
becomes 

tl ^ t^^m^ hm n^2^ 

° ° AO^O ^(f]2 _^2)2 _^^2^2f|2' 
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In the following, we suppress the labels in the frequencies, since they are continuous quantities. 
We start by defining a function W{z), 



oc 2 2 

W{z) ^ - u;^ + Y, '^^ 



We find that the fi's are the roots of W{z). Using 77^ = 2gAuj, we have in the continuum limit. 



duj 



W{z) = z^ -Cj'^ + 2gz' 



For complex values of z the above integral is well defined and is evaluated by using Cauchy theorem, to be 



W{z) 



z^ + ig-KZ — w^, Im(z) > 
z^ — igirz — ili^, Im < . 



We now compute /oo(i) — X]^^o(^o)^^ * which, in the continuum limit, is given by 



/oo(i) 



We find that. 



^ °^ " Win) ' 

and since the fl's are the roots of W{z), we write Eq. ([55]) as 



/oo(t) 



(85) 



(86) 



(87) 



(88) 



(89) 



(90) 



where C is a counterclockwise contour in the z-plane that encircles the real positive roots of W{z). Choosing a contour 
infinitesimally close to the positive real axis, that is z = a — ie below it and z = a + ie above it with a > and 
e ^ 0"'', we obtain 



Mt) = - 



daae 



1 



1 



W{a-ie) W{a + ie) 



In the limit e ^ 0+, Eq. ([ST]) gives W{a ± ie) = — ± igira which leads to 

fooit) = Ci{t;u),g) + iSi{t;ui, g), 

where 

cos(at) 
(a^ ~ tj2)2 + TT^g^a'^ ' 
o? sin(Q;t) 



G\{t\^,g) = 25 / da 



S\{t\uj,g) = -2g 



da 



(a^ — LJ^)^ + TT^g^Q!^ 



(91) 
(92) 

(93) 
(94) 



Notice that Ci{t = 0;u!, g) = 1 and Si{t = 0;ii>,g) = 0, so that foo{t = 0) = 1 as expected from the orthonormality 
of the matrix (tp. The real part of foo{t) is calculated using the residue theorem. For = tD^ — n^g'^/A > 0, which 
includes the weak coupling regime, one finds 

Ci{t;LU,g)^e-^s'/^ 



cos(Kt) — — ^ sin(Kf) (k > 0). 
2k J 



(95) 



Although Si{t;u},g) cannot be analytically evaluated for all t, however for long times, i.e. 1/^^, we have 

S^{t;u,g)^^ (i»l). (96) 
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Thus, we get for large t 



l/ooWl 



2 ^ ^-.gt 



cos[Kt) — — sin(Ki) 
2k 



+ 



(97) 



Next we compute the quantity fak{t) — X^J^o ^o^I-^ * ^'^^ continuum Umit. It is 



where 7y = ^y2gALu. Taking the same contour as that used to calculate fooit), we obtain 



da 



W{a - ie)[{a ~ itf ~ ^^1 W(a + it){(a + ief ~ uj'^ 



Thus, taking e ^ 0+ fou{t) is written as 

foLo{t) = wVAw [C2(ti^,t;(:),g) + iS'2(w,i;tl',5)] 

where 

cos(at) 



C2(t^,t;w,g) - (25)5 
S2{uj,t;uj,g) = -{2g) 



da- 



o? sin(Q;t) 



da. 



(99) 

(100) 

(101) 
(102) 



Notice that the integrals defining the functions C2 and S2 are actually Cauchy principal values. 
The function C2 is calculated analytically using Cauchy theorem; we find 



C2{uj,t;uj,g) 



2ff 
wg 



-Tigt/2 



2k (w^ — a;2)2 + TT^g'^uj^ 
■ngu 



(w^ — a;2)2 + TT^g^CJ^ 

sin 



cos nt 



sin tjt 



(103) 



(cj^ — Lo'^Y + Tr'^g^uj^ 

The function ^2 cannot be evaluated analytically for all t, it has to be calculated numerically. For long times, we have 



(104) 



In the continuum limit, we get the average of the particle ocupation number, 

n'o{t) = [Cf{t;iu,g) + Sf{t;cu,g)]n'o+ / dcocj^ [C^{u;,t;iu,g) + Si{iu,t;ii,,g)]n'{u;) , (105) 



where n'{uj) — l/{e^'^ — 1) is the density of occupation of the environment modes, the functions Ci and C2 are given 
by Eqs. ([95)1 and (|103p while the functions Si and S2 are given by the integrals Eqs. (|94p and (|102p . respectively. In 
Fig. [2] we display the behavior in time for no = 1, a; = 1, /3 = 2 and g = 0.1; {t > 1). 

The important point, that is seen from Fig. [T] and Fig. [2] is that, for long times, both the bare and dressed 
ocupation numbers of the particle approach smoothly to asymptotic values which are ~ 0.160. Moreover these values 
are expected on physical grounds, being slightly higher than the one obtained from the Bose distribution at the 
equilibrium temperature of the reservoir. In fact, taking (3 — 2 and cj = 1, as used in the plots, one has, 

noo(<I') = l/(e'^" - 1) = 0.156 
Therefore both methods, and in particular our dressed state formalism describes correctly the thermalization process. 
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FIG. 2: Time behavior for no(t) given by Eq. (fTOSl) . for {t > 1), no ^ 1, lj ^ 1, l3 = 2 and g = 0.1 



VI. FINAL REMARKS 



We have considered a linearized version of a particle-environment system and we have carried out a non-perturbative 
treatment of the thermalization process. We have adopted the point of view of renouncing to an approach very close to 
the real behavior of a nonlinear system, to study instead a linear model. As a counterpart, an exact solution has been 
possible. This realises a good compromise between physical reality and mathematical reliability. We have presented 
an ohmic quantum system consisting of a particle, in the larger sense of a material body, an atom or a Brownian 
particle coupled to an environment modelled by non-interacting oscillators. We have used the formalism of dressed 
states to perform a non-perturbative study of the time evolution of the system, contained in a cavity or in free space. 
Distinctly to what happens in the bare coordinate approach, in the dressed coordinate approach no renormalization 
procedure is needed. Our renormalized coordinates contain in themselves the renormalization aspects. As far as the 
thermalization process is concerned from a physical viewpoint, both bare and dressed approaches are in agreement 
with what we expect for this process. For long times, all the information about the particle occupation numbers 
depends only on the environment. Both curves in Fig. [1] and Fig. [5] approach steadly to an asymptotic value of the 
bare and dressed ocupation numbers of the particle, which is the physically expected one at the given temperature. 
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